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Abstract 

The mapping of steady-state nonequilibrium dynamical mean-field theory from the lattice to the impurity is described in 
detail. Our focus is on the case with current flow under a constant dc electric field of arbitrary magnitude. In addition to 
formulating the problem via path integrals and functional derivatives, we also describe the distribution function dependence 
of the retarded and advanced Green's functions. Our formal developments are exact for the Falicov-Kimball model. We also 
show how these formal developments are modified for more complicated models (like the Hubbard model). 
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1. Introduction 

Dynamical mean-field theory (DMFT) began in 
1989, wlien Metzner and Vollliardt suggested the 
large-dimensional limit (with an appropriate rescaling 
of the hopping integral) as a simplifying limit for the 
many-body problem that nevertheless still included 
the important competition between minimizing the 
kinetic and potential energies of the system [1]. These 
ideas were implemented into a complete DMFT by 
Brandt and Mielsch with the exact solution of the 
Falicov-Kimball model at half filling [2]. Since that 
time, nearly all equilibrium many-body models have 
been solved in the large-dimensional limit with DMFT. 
Recently, progress has been made on nonequilibrium 
extensions of DMFT [3,4,5,6,7] for the Falicov-Kimball 
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model and for the Hubbard model [8]. In situations 
where the perturbing fields (or the time-dependent 
part of the Hamiltoniari) maintains the translational 
invariance of the lattice, the self-energy remains uni- 
form and local, as follows from the Langreth rules [9] 
and the perturbative expansion in the hopping by 
Metzner [10]. 

In this work, we examine the steady-state limit of 
nonequilibrium DMFT. We start the system in equi- 
librium at an inverse temperature /3 and wait a long 
time after the field or time-dependence has been turned 
on, so that the system has reorganized itself to the 
long-time response of the driving fields (or time depen- 
dence) [11,12]. We are inherently assuming that the 
long-time limit of the Hamiltonian is different from the 
original equilibrium Hamiltonian. It need not have any 
explicit time dependence in this limit though. We will 
show that the Keldysh formulation of the many-body 
theory {ignoring the third (imaginary) branch of the 
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contour described by Wagner [13]} produces the exact 

steady-state solution for the Fahcov-Kimball model; we 
will also discuss modifications for the Hubbard model 
(and more complicated models). Our approach works 
with functional integrals and derivatives. Most of the 
techniques are completely straightforward, but these 
details have not appeared in the literature yet, and pro- 
vide useful insights into the solution of nonequilibrium 
problems via DMFT. 



2. Formalism for the Falicov-Kimball model 

We will focus our efforts first on the Falicov-Kimball 

model, whoso Hamiltonian [15] is the following [in a 
uniform electric field described by a uniform vector 
potential E = -dA{t)/cdi\: 

Hfk = ^{e[k - eA{t)/hc] - njc^c^ + ^ X/ ^1^*^'' 

k i 

(1) 

where we employed the Peierls' substitution [14]. 
Here, and cj^ destroy and create a spinless itinerant 
fermion with momentum k, e(k) = —limd^ca(t*/y/d) 
Ylt^i ^'^^ ^® bandstructure on an infinite- 

dimensional hypercubic lattice, fj, is the itinerant- 
electron chemical potential, U is the conduction- 
electron-local-electron interaction, Wi is the number 
operator for the spinless localized electrons at site 
i, and we use the real-space basis for the itinerant 
electron operators in the second term of H. In the 
following we set h = c = t* = 1. 

We are interested in the steady-state limit of the re- 
sponse, which can be determined in the limit where the 
vector potential is turned on in the infinite past (but 
after the system has fully equilibrated into an equilib- 
rium distribution with inverse temperature iu other 
words, we choose A{t) = limto_>-oo 0{t—to){t—to)E. In 
nonequilibrium physics there are two different Green's 
functions that are required to describe the system. One 
is the retarded Green's function and the other is the 
Keldysh Green's function. The former describes how 
the quantum density of states (DOS) varies with en- 
ergy, while the latter describes how the electrons axe 
distributed amongst those states. Naively, one would 
assume that the quantum states would be indepen- 
dent of the initial temperature, although they will de- 
pend on the electric field strength, the fillings of the 



particles, and the interaction energy. But it is well 

known in equilibrium many-body physics that inter- 
acting systems often display some temperature depen- 
dence to the DOS. Noninteracting problems, and the 
Falicov-Kimball model are exceptions, however, where 
the DOS is independent of how the electrons are dis- 
tributed amongst the states. It is the goal of this work 
to elaborate on this situation in nonequilibrium cases. 

Our starting point is to define the contour-ordered 
Green's function on the two-branch Keldysh contour, 
which runs from t — — oo to f = cx) and back. We define 
it only for the case of an impurity in a time-dependent 
dynamical mean field denoted Ac: 

Gc{t,t') = -iTr,f{TcpimpS{Xc)c{t)cHt')}/Z, (2) 

with pimp the steady-state density matrix for the impu- 
rity, the operators denote time evolution in the Heisen- 

berg picture with respect to the impurity Hamiltonian 
Himp = —ne'e + UwiC'c, the S-matrix (evolution op- 
erator) satisfies 



i ( dt ( dt'c\t)\c{t,t')c{t') 
Jc Jc 



S{Xc) = Tcexp 

.In .In 

(3) 

2 = Tr c f pimpTc S (Xc) is the partition function, the 
time ordering is with respect to the ordering along the 
contour, and the trace is over the four states of the 
spinless conduction and localized electrons. The long- 
time-limit density matrix is a priori unknown (while 
the lambda field is determined self-consistently with 
the DMFT algorithm). In practice, one needs to evolve 
the full many-body system from the time the field is 
turned on until the long-time limit is reached (and use 
the so-called restart theorem [16] to determine it). We 
will see, however, that the density matrix does not en- 
ter into the calculation of retarded or advanced quanti- 
ties for the Falicov-Kimball model, so we defer further 
discussion at the moment. We will work with Wigner 
coordinates [17] of average T = {t + t')/2 and relative 
trei = t — t' time below. 

We begin our discussion by considering a problem 
that has no localized particles, so we do not take the 
trace over the /-electrons, but instead, simply set Wi = 
0. The full solution, including the trace over /-particles 
can be easily constructed from this solution, as we show 
below. It is convenient to break up the two-branch 
contour into a + branch, where the time increases 
from — oo to -|-oo and a minus branch, where time de- 
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creases from +00 to —00. Then the evolution operator 
in Eq. (3) can be rewritten as 



S{Xc)=Tcexp -i dt dt'l (4) 

J — 00 J — oc I 

- cl {t)X> {t, ?)c+(?) + cL {t)X^{t, F)c_ (?) 



where the fermionic operators with the + or — sub- 
script hve on the corresponding time branch; the mi- 
nus sign enters from the change of direction in how 
the time evolves on the two different branches. When 
both time arguments lie on the -|- branch or the — 
branch, we have time-ordered (T) or anti-time-ordered 
(T) objects, respectively. When one is on the -|- and 
one on the — , we have have the lesser (<) or greater (> 
) functions. To be concrete, the four different Green's 
functions are defined as follows (we suppress the trace 
over the /-electrons, which is needed for the Falicov- 
Kimball model, but is neglected for the moment): 

G^{t,t') = -iTrc{TtpimpS{Xc)c+{t)cl{t')}/Z, (5) 

G^{t,t') = -iTrc{TtPimpS{Xc)c_{t)cUt')}/Z, (6) 

G<{t,t')=iTrc{pimpS{Xc)cl{t')c_{t)}/Z, (7) 

G>it,t') = -iTr,{pir„^SiXc)c+{t)ciit')}/Z, (8) 

where the t and t subscripts denote time-ordering or 
anti-time-ordering, respectively. The retarded and ad- 
vanced Green's functions are defined to be G'" = G^ — 
G^ and G" = —G^ + G^, which can be written as 

G"'(t,t') = -ie{t - t')Tr, [pim.pSiXc){c{t), ct (?)}+] 

(9) 

G^{t,t') = i6{t' -t)Ti, [pimrS{Xc){c{t),c^{t')}+] 

(10) 

in terms of the operators. Similarly, the Keldysh and 
anti-Keldysli Green's functions are defined to be G^ = 
G> + G< and = -G^ - G'^ + G> + G<, or 

G^(t, t') = -iTrc {pimpS{XcMt), ct (t')]- } , (11) 
G^(t, t') = -iTrc {pimpS{Xc){[c{t), ct (t')]- 

-[c(t),ct(t')]-}}. (12) 



Note that the anti-Keldysh Green's function vanishes, 

but we need its functional form in order to take func- 
tional derivatives for different Green's functions below. 
Using the definition of the evolution operator, we im- 
mediately find that the Green's functions can be found 
as functional derivatives of the partition function. The 
explicit relations are 



5X<{t',t)' 



for the time-ordered, anti-time-ordered, lesser, and 

greater Green's functions, respectively. We will dis- 
cuss the alternative retarded, advanced, Keldysh and 
anti-Keldysh basis below. Now that we have all of 
these definitions, we can actually solve explicitly for 
these Green's functions using the equations of motion 
(EOMs). We take the defintion of the Green's function, 
and differentiate with respect to time; where appropri- 
ate, one needs to take into account the time ordering, 
which brings down terms proportional to the A fields. 
This procedure is straightforward to complete, and 
the end result is the following: 

idtG'^it, t') = 5{t - t') - p.G'^it, t') 

/oo 
dtX^ (t,i)G'^ (i,t') 
-00 

/OO 
dtX<{t,i)G>{t,t'), (15) 
-00 

- idtG^it, t') = S{t - t') + ixG^it, t') 

/OO 
dtX>{t,t)G<{t,t') 
- 00 

/OO _ 
dtX^{t,i)G'^{i,t'), (16) 
-00 

/OO 
dtX^{t,i)G<{i,t') 
-00 

/OO 
dtX'^{t,t)G^{t,t'), (17) 
-00 

/OC 
diX^{t,t)G'^{i,t') 
-00 

/OO 
diX^{t,t)G>{i,t'), (18) 
-00 
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with similar equations for derivatives with respect to 
t' which we do not write down explicitly. 

On the lattice, one can exactly solve for the anal- 
ogous contour-ordered Green's functions in the case 
where ?7 = [18]. One can sec from the exact solution, 
that the noninteracting retarded Green's function be- 
comes average time independent for long times after 
the field has been turned on. We make the ansatz that 
the retarded self-energy is also independent of average 
time in the long-time limit (which is consistent with 
gauge- invariancc arguments [19,20,7], but is an inde- 
pendent ansatz). Then, it is straightforward to show 
that the interacting retarded Green's function is in- 
dependent of average time. In general, we can only 
show that the average time dependence of the contour- 
ordered Green's function is periodic in average time 
with a period given by the Bloch period 2ii / E [8]. This 
allows us to make a continuous Fourier transform with 
respect to the relative time (frequency dependence a;) 
and a discrete Fourier transform with respect to the 
average time (Fourier components tilub ~ nE) for all 
of the different components of the contour-ordered ob- 
jects (time-ordered, anti-time-ordered, lesser, greater, 
and Keldysh); the retarded and advanced components 
depend only on the continuous frequency u). 

In order to find the Green's functions, we need to em- 
ploy the simplification from the retarded and advanced 
Green's functions, which depend only on the relative 
time, or frequency oj. Hence, we need the EOM for the 
retarded and advanced functions, which can easily be 
shown to satisfy 

/oo 
diy{t-t)G''(t), (19) 
-oo 

/oo 
dtX'{t-t)G''{t), (20) 
-OO 



which are solved by 



uj + ^ — 



UJ + H - ' 
(21) 

where we have solved the problem in frequency space 
after a Fourier transformation. Next, we can show, by 
using the identities that relate the different quantities 
(T, T, <, >, r, and a) to each other, that the EOM for 
becomes 



/oo 
diX'{t-^G^{t,t') 
-OO 

/oo 
dtX<{t,t)G''(t-t'), (22) 
-OO 

since = A'' -|- A< and G"^ = G< -G". Substituting 
in the appropriate Fourier expansions then yields 



(u! + ^aJB)G^(nuJB,u!) = \^{uJ + —wb)G'^{iuob,uj) 



+ \'^(nujB,co)G'^{uJ — 2^3)- 
This can be directly solved to yield 



G^{nwB,oj) = 



D{UJ+%UJB,U>- ^uib)' 



(23) 



(24) 



with 



71 Tl Tl Tl 

D{uj + -u>B,0J - -^'^b) = [''^ -I- -^^B + 1^- y{uj + -2^b) 



Th Tl 

X [w - —ljb + M - A''(a; - -ws)]. 



(25) 



Similarly, we have 

G^{nujB,i^) = 



D{UJ+%U>B,UJ- ^ujb)' 



(26) 



and G"^ {tiujb,(^) = G'^{nujB,(^) + SrMG^{uj) and 
G'^{nuiB,u}) = G^{nu3B,u!) — 6„oG°-{lu); the Kro- 
nnecker delta functions enter because the retarded 
and advanced Green's functions are independent of 
average time. One can directly show that 

X'"(nuJb,uj) = X^inuJbjUj) = A^(nu;6,w) = A^(nc<;i,,c<;), 

(27) 

for n 7^ 0. Now replacing the retarded and advanced 
quantities in the denominator by the time-ordered, 
anti-time-ordered, lesser and greater quantities, gives 

D=[LU+!^OJB+fi-X^(0,LU+'^iJB)] (28) 

X [uj - -ujb + ^^ - X (0,ijj - -ljbJI 

-h A<(0,w + ^ujB)[-nojB + X>iO,u) - ^ojb)], 

which involves just the zeroth component of the Fourier 
series terms for the dynamical mean fields. 

Finally, we can now solve the functional differ- 
ential equations, which take the form G{nuJi,-u>) = 
±5\n2/6X{—moB,'jj), and determine the partition 
function as 
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lnZ = J duj In ^{oj + - {Q,u})){ui + 

A^(0,t<;)) + A<(0,cj)A>(0,a;)j +C 



Using these solutions, we can integrate them to find the 

partition function, but we need to introduce appropri- 
(29^tc anti-Keldysh fields to give us the final functional 
form, because we set that field to zero in all functional 



D{u>+ fa)s,w - fws) 



where we fix an undetermined constant C by equating 
with the noninteracting result. It is a straightforward 
exercise to show that the functional derivatives of this 
partition function yield the Green's functions. 

Our next step is to repeat these results for the alter- 
native r, a, K, K basis. To start, the evolution operar- 
tor becomes 



S{\c) = Tt exp 



I ^ >/ — oo J — c 



dt' 



[cl{t)-cl{tW{t,t')\c+{t')+c.{t')] 
+ [c\{t)+cl{tW{t,t')[c+{t') - C-{t')] 
+ [cl{t) - cl{t)]\''{t,t')[c+{t') - C-{t')] 



where we have added a now field A which will bo set 
equal to zero for all physical matrix elements that we 
evaluate. Using the definitions for the retarded, ad- 
vanced, Koldysh and anti-Kcldysh Green's functions 
shows that we can determine them via functional 
derivatives 

where all derivatives must be evaluated with A^ = 
after taking the derivatives. Since we have already de- 
termined the retarded and advanced Green's functions 
when we solved for the Green's functions in the other 
basis, we need only solve for the Keldysh and anti- 
Keldysh Green's functions, which can be easily deter- 
mined (via the EOM, or using the relation between the 
lesser and greater Green's functions) to give 



G (nu)B,ijj) = 



\'^{noJB,uj) 
D{lj+^ujb,uj-^ujb)) 



(33) 
0. (34) 



derivatives we evaluate. The end result is 



In 



Z = j dw\n[{uJ + |J,-W{uJ)){uJ-\-^l-\'^{uJ)) 

- A^(0,a;)A^(0,c^)] -hC 
X - A^(nws, w)A*^(— nwB, w) 



(35) 



This completes the derivations for the two equivalent 
bases for the impurity model with no localized elec- 
trons. 

We now generalize for the case of the Falicov-Kimball 
model, where we include the trace over the localized 
particles in all of the relevant expectation values. The 
changes are straightforward to work out, and we re- 
port them only in the r, a, K, and K basis. Define the 
effective partition function by (noting that A'' = A"*) 



Zo{n) = C exp 



/ 



dui 



Iw-FAt- A'-(w)|2 



(36) 



where C = Tr cPimp, is the noninteracting partition 
function for the steady state in the absence of the time- 
dependent fields. Then the full partition function is 



(37) 



with Ef the localized electron Fermi level, adjusted to 
give the correct average filling of the localized particles. 
The average filling satisfies 



1 — Wo 



e~'"'fZo{n-U)_ 



Wo 



2^0 in) 



■ (38) 



The Green's functions take the same form as they had 

for the case with no localized particles, except we now 
have the sum of two terms: one, weighted by wo, which 
is the result we have found above, and one, weighted 
by Wl, which has the same form as the result we de- 
rived above but with n —> n — U. The retarded Green's 
function takes the same functional form as it has in 
equilibrium, namely 



G'-(o;) 



■ Wl 



Wl 



(39) 



The Keldysh Green's function is more complicated: 
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{nu)B,ijj) = X'^ (nujB,oj) 



+ 



1 — Wl 



D{u3 + |wb,w - fWs)U^;u-f/ 



(40) 



Note that the retarded Green's function appears to 
be independent of the Keldysh Green's function, and 
hence has no exphcit dependence on the distribution 
function, and thereby, no dependence on the steady- 
state density matrix, but one needs to carefully check 
to see whether the functional derivative of wi has any 
dependence on A'^^. Indeed, it is very easy to show that 
there is no such dependence, because any functional 
derivative of the partition function with respect to 
will bring down a factor of A^^, which is set equal to 
zero, so we immediately learn that 



SG^t,t') 1 



(43) 



x{ - 2c{t)Tt[c^t')c^{t')c{t)] + 2Ti[c\t')c^i')cit)]cit) 
+Tt[c^ {t')cit)]Ti[c^t')c{t)] + %[c\t')c\t')]T,[cit)c{t)] 
+%Ht)c' {t'mic' {t')cit)] - % [ct (i')c(t)]r,[c(t)ct 
-Tt [c{t)c{i)]Tt[c^ (t>t (?)] - Tt [ct (t')c(f)]Tf [ct {t')cit)] } . 

We know that this correlation function vanishes for the 
Falicov-Kimball model, but it is nontrivial to show that 
this is so, and it docs not appear obvious at all from 
the operator expectation value above. 

Note that we cannot make any further progress on 
evaluating the Keldysh Green's function in the steady 
state without knowing what the density matrix is, and 
we do not know this explicitly. 



<5'"G"(a;) 



(5A^(niWB,wi) . . .6X'^{nm.U}B,U>',n) 



= 0, (41) 



for all powers m. This is the generalization of the proof 
that the conduction-electron density of states for the 
Falicov-Kimball model is independent of temperature in 
equilibrium to a proof of the distribution-function inde- 
pendence of the conduction- electron density of states in 
the nonequilibrium "steady state". Of course, this proof 
only holds in the case of a canonical distribution of the 
heavy particles, where their filling is kept constant as 
the temperature or the nonequilibrium driving fields 
are varied. 

In order to understand this result more fully, wc 
examine explicitly the functional derivative of the re- 
tarded Green's function with respect to the Keldysh 
dynamical mean field in the time basis. We find that 



SG''{t,t') _ 1 
<5A^(f',t) ~ 4 



T^cfTcPimpS{Xc) 



x[clit')-cUmC+{t)+C.it)] 

X [CU?) - Cl(?)][c+(t) - C.{t)]/ZFK. 



(42) 



The -I- and — subscripts denote which branch of the 
Keldysh contour the time arguments lie on. This ex- 
pression can be simplified somewhat to the form 



3. Formalism for the Hubbard model 

In the case of the Hubbard model, both particles 
can move, so we have an evolution operator given by 
the product S{Xj-;)S{Xj-,) in all of the expectation val- 
ues that we need to evaluate; the | electrons are the 
old c electrons and the J. electrons are the old / elec- 
trons. Now we can set up the formalism for calculating 
the Green's functions directly, just as for the Falicov- 
Kimball model, but unfortunately, in this case, the 
equations of motion cannot be solved explicitly, so we 
cannot determine an explicit formula for the partition 
function anymore. These expectation values can be cal- 
culated with the numerical renormalization group, if 
wc have the explicit operator form for the density ma- 
trix, but as we discussed above, this is not the case. 

What we do have is an explicit operator average that 
tells us the dependence of the Green's function on the 
distribution function (via the Keldysh dynamical mean 
field). It is just the result in Eq. (43), generalized to 
include spin indices a on the operators with unbarred 
time arguments, and a' on operators with barred time 
arguments. If that correlation function could be eval- 
uated for an (approximate) solution of the retarded 
Green's function for the Hubbard model [8], then we 
could determine how strong the distribution function 
dependence of the retarded Green's function was for 
that model. Unfortunately, this does not appear to be 
a simple task at this point in time. 
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4. Conclusions 



[7] N. Tsuji, T. Oka, H. Aoki, Phys. Rev. B 78 (2008) 
235124. 



In this work, we have shown a number of exphcit 
details for the impurity problem in nonequilibrium dy- 
namical moan-fiold theory when one approaches the 
steady-state response. For the Falicov-Kimball model, 
we can advance the functional analysis to the point 
where we can prove that the retarded (and advanced) 
Green's function does not depend on the distribution 
function. In addition, we derived a two-particle corre- 
lation function, that would need to be evaluated for 
the general case, which determines the strength of the 
distribution-function dependence of the retarded (and 
advanced) Green's function. While wc can explicitly 
show this correlation function vanishes for the Falicov- 
Kimball model, we do not know how one would eval- 
uate it for the Hubbard model, particularly because 
we do not know what the steady-state density matrix 
is for the model. Nevertheless, we anticipate that it 
is not large, as wo expect the steady-state nonequilib- 
rium density of states to depend only weakly on the 
distribution function. 
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